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For the last decade the reaction-diffusion system A + B ^ 0, where unhke species A 
and B diffuse and irreversibly react in the bulk of a (i-dimensional substrate, has acquired 
the status of one of the most popular objects in nonequilibrium statistical physics [1]. 
Two problems have been investigated most intensively: (i) dynamics of the A + B ^ 
annihilation in an infinite system with initially homogeneously (randomly) and equimo- 
larly distributed reactants, in which below the critical dimension dc = 4: a dynamical 
clustering develops (fluctuation-induced like-particle domains formation) and, as a result, 
an anomalous reaction deceleration arises (Ovchinnikov and Zeldovich, 1978 [2], Toussaint 
and Wilczek, 1983 [3]); (ii) behaviour of the dynamic reaction front in an infinite system 
with initially spatially separated reactants (Galfi and Racz, 1988 [4]), and structure of 
the steady state front in a finite system, at the ends of which are injected equal currents 
of A's and B's particles (Ben-Naim and Redner, 1992 [5], Cornell and Droz, 1993 [6]). 

In 1994 it was discovered (Shipilevsky, 1994 [7]) that in another wide class of RD 
systems, where reaction and diffusion are spatially separated (i.e. reaction proceeds on the 
surface of the medium and diffusion proceeds in its bulk) the interplay between reaction 
and diffusion acquires qualitatively new features and leads to the threshold self-organizing 
dynamics of the A + B ^ 0. It has been found that once particles A and B diffuse at 
different mobilities from the bulk of finite medium onto the surface and die on it by the 
reaction A + B — >■ 0, there should exist some threshold difference in the initial numbers of 
A and B particles, A^, above which the loop of positive feedback is "switched on" and the 
process of their death, instead of usual deceleration, starts to accelerate autocatalytically. 
Recently, I have shown [8] that deceleration-acceleration transition is only a prelude to 
much more nontrivial dynamical effects, which take place far beyond the acceleration 
threshold. It has been demonstrated that in the diffusion-controlled limit A — oo a new 
critical phenomenon develops, which arises as a result of self-organizing explosive growth 
(drop) of the surface concentrations of, respectively, slow and fast particles [concentration 
explosion) and manifests itself in the form of an abrupt singular jump in the desorption 
flux relaxation rate. The aim of this report is to give an insight into the key features of 
this remarkable phenomenon which I call annihilation catastrophe. 

Model. I consider a model, in which species A and B are supposed to be initially 
uniformly distributed in the bulk of infinitely extended slab of thickness 2i. Both species 
diffuse to the surface (X = ±i) and desorb as a result of surface reaction Aads + Bads 
AB — > 0. The rate of reaction is proportional to the product of surface concentrations 
/ = kpAPB = t^CAsCBs, where the rates of surface-subsurface exchange are supposed 
to be rather high, so that quasiequilibrium pi = fiCis is always sustained. The system is 
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effectively one dimensional and the boundary conditions are determined from the equality 
of diffusion and desorption / flux densities at the surface \s= /, i.e. it is assumed 
that the surface layer capacity can be neglected. Thus, introducing index (/ieaf?/) 
for slower diffusing species and index " L" {light) for a faster one, the problem reads ( by 
symmetry I consider the interval [0, £] only) 

dh/dr = V'/i , dl/dr = {l/p)VH, (1) 

Vh 1= (l/p)V/ |,= -Kh, (2) 

with V(/i,/) |x=o= and the initial conditions h{x^{)) = ho and l{x,0) = Iq. Here 
h{x, t) = ch/c^ and /(x, r) = cl/c* are the reduced concentrations, V = d/dx, x = X/i e 
[0, 1] is the dimensionless coordinate, r = Dnt/i"^ is the dimensionless time, p = Dh/Dl < 
1 is the ratio of species diffusivities, and c* = Dh/k£ is the characteristic concentration 
scale, from which diffusion starts to play an essential role. Boundary conditions (2) 
reflect the fact that particles disappear in pairs only, i.e. J = hsh = — < h >= — < I >, 
where J = ///* is the reduced desorption flux density and /* = kcI is its characteristic 
scale, therefore < h > — < I >— A = const., i.e., the excess amount stays "inert" in 
the bulk ( here < h >= Jo hdx = J^h/J^* and < I >= Jo Idx = Nl/.N'* are the total 
reduced numbers of particles in the bulk per unit of surface and A/"* = = Dh /k is the 
characteristic particle number scale). I consider here the annihilation dynamics for A > 
when the pair number N is dictated by the L particles number N{t) = Mpair/^* = ^l^t) 
so that in the final state all L particles disappear N{oo) — l{x, oo) = and H particles 
are distributed uniformly with concentration h{x, oo) = A. I am interested first and 
foremost in the diffusion-controlled limit A^o ^ oo (note that diffusion-controlled limit 
No = Iq = cl(0)/c* — >• oo means = h£/Dh — >■ oo). Moreover, I focus on the 
limit no = A'o/A ^ 1 when the initial numbers of H and L particles are almost equal, 
r = ho/lo = {1 + no)/no ~ 1. Precisely in this limit the particles interplay is most 
"intensive" and gives rise to the explosive dynamics, that is my prime interest here. 

Transition point and self-accelerating annihilation. Assuming A^^q ^ oo and 
TZ = Tyjp < 1, I show [10], that for no ^ 1 and 7 = 1 — 7^ l/i/iVo at times 
T ^ Th — 1/Nq ^ the annihilation is followed by a sharp drop of hg, crossing over to the 
if-diffusion-controlled regime, where the desorption flux density is governed by the rate 
of H particles diffusion to the system's surface. In the range <C r ^ 1, when diffusion 
of H particles proceeds in a semi-infinite medium regime, J ~ hoj s/tvt ~ Jq/ ^/T^iTh/T)^^"^ 
and for surface concentrations I find 

hs = (r/7V^)(l + ^ + • • •), ls = /o7(l - ^ + • • •), 

where q — 2r^T/'K at p ^ 1. At r > 1 the slow power-law stage crosses over to the 
exponential relaxation stage J = Ae'''^°'^{l -\- e~^'^°'^ -|- • • •) with amplitude A — 2ho, 
and I show that in the limit ^ — > 00 the surface concentrations achieve universal (A- 

independent) asymptotics 

= A,(l + + e-^°^ + ■■■), {A/A,)e-'^°-{l -</> + •••) (3) 
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where |0| decays as 0{e *^'^°^) in the range p < Pc = 4/9 and as 0{e ^'^^/^ i^o)t^ ^j^g 
range Pc < p < 1 [10], and Ac obeys the expression [7,9] 

Ac = ^Juj^isj[i{^Juj^) 

{(jJq = 7r^/4 being the main eigenfrequency of the diffusion field relaxation under absorbing 
boundary conditions) . This important result implies that in the diffusion-controlled limit 

A ^ oo the critical asymptotics hg = hg = Ac appears at any A as a precursor of long- 
time asymptotics h^^ \^^^= A, and precisely this intermediate asymptotics selects 
the critical value A = Ac = /i^ below and above of which the relaxation dynamics of 
hg is qualitatively different: at A < Ac the surface concentration, hg, escaping from the 
asymptotics h^, decreases, reaching its limiting value A from above, whereas at A > Ac the 
surface concentration, hg, escaping from the asymptotics h^, increases, reaching its limiting 
value A from below [7]. With neglect of rapidly decaying terms of 0{e~^'^°'^,e^'^ /p-^o)t^ 
the exact asymptotic expansion in power of 1/^4 has the form [10] 

hs = Ac(l + A + • ■ •), ls = {A/A,)e-^°^{1 - A + ■ • ■), N= (^/cc;o)e-"°^(l - /.A + ■ ■ ■), 

J = ^e-'^o"(l - v^tanh ^/^{A - Ac)"^^ + •••), (4) 

with n = u!o/Ac ^ 1—p and A = Ac(A — Ac)v4.~^e'^"^. These expressions and more detailed 
analysis in Ref. [7,9] suggest that at the crossing of the critical point A = Ac not only the 
sign of hg changes, but a growth of hg (and with it the whole annihilation process) begins 
to accelerate autocatalytically, kg — cuo{hg — Ac) + • • •, as a result of " switching-on" of 
the loop of positive feedback: diffusion- induced growth of hg accelerates the drop of Ig, 
the drop of Ig accelerates the growth of hg and so on {self-accelerating pairs death). This 
phenomenon, interpreted in [7] as a new type of self-organization, develops the brighter 
the stronger A exceeds Ac, and in the hmit A ^ oo it results in the explosive dynamics 
- annihilation catastrophe [8]. 

Annihilation catastrophe. Scaling and universality. According to (4) in the 
limit A — s> cxD at the self-acceleration stage the rates of growth Qhs = +d\nhg/ dr and 
relaxation VLls = —dlnlg/dr of surface concentrations are unambiguously defined by the 
current value of the reduced pair number ^(t) = N{t)/A and are mutually "compen- 
sated" so that with an accuracy to vanishingly small terms 0(l/n^A) the dcsorption 
flux relaxation rate rj^ = —d\n J/dr = n^g — ^Ihs = c<Jo + ■ ■ ■ is sustained constant. The 
central result of [8] is my finding, that at A — (X) the relaxation rate rj^ is sustained 
constant up to the time moment of reaching some critical pair number n^^ip), in whose 
vicinity 5t oc A~^/^ — > the tJ^ abruptly jumps from tJ^ — ujq to tJ^ — > oo. I consider 
here the two consecutive stages of this phenomenon. 

i) Concentration explosion. Following the work [8], 1 first show that condition tJ^ = luq 
in the limit A ^ oo leads automatically to the arising of a finite-time singularity. Taking 
J = Ae~^°'^ and satisfying (1),(2) to the leading order in A, I find 

j = J/A = -n = uo{l + n), kg = Ac{l + n)/{n - n^), = /xA(n - n*), (5) 

where = (1 — /x)//i ~ p/(l — p). It is seen that kg diverges as n — >• and at the 
point of singularity exactly = Ac. Prom equality J* = 2hoe~'^°'^* — AAc in accord with 
the numerical data of [10] we obtain 
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r, = (4/7r2)ln[2(l + no)/Ac]. 



Taking the origin of time at the point r^, i.e. introducing a relative time T = t — t^, 
from (5) we come to remarkably simple laws n = (l/^)e'^"l^l ~ 1-J = A^.e'^"'^', /i^ = 
AJ{1 - e-ol^l),/, = A(e-ol^l - 1), and, hence, Qhs = ^i{hs - A^) = u;o/(e"»l^l - 1), 
whence it follows that for any p < 1 at the "distance" from the singularity point |T| ~ 
the dynamics of growth of hs changes drastically: the initial exponential increase in the 
growth rate VLhs — <^oe~'^°''^'(f^Ls ~ <^o) turns at ujq\T\ -C 1 to a singular one 

^Hs^^Ls^ IJ'hs^l/lTl |r|^0, (6) 

i.e., concentration explosion develops in the system. The striking feature of this phe- 
nomenon is that at times ujq\T\ -C 1 the desorption flux practically reaches its criti- 
cal value J*, therefore the concentration explosion develops, in fact, at a "frozen" flux 
hsls = J*,hs = l/fi\T\,ls = fJ'J*\'^\j the consequence of which is synchronization of the 
growth and relaxation rates Qhs = ^Ls = during the explosion. A singular growth of 
hs at the background of the diffusion flux J*, "frozen" on the explosion time scale, gives 
rise to a singular growth of the diffusion " antiflux" , the calculation of which yields [8] 
~ — 1/A|T|^/^. As a result, in the desorption flux relaxation rate rj^ — u;o + [rj^] 
the singularly growing term arises 

[tJ'] ~ l/A\Tf/', (7) 

whence it follows that the scenario described above holds down to \T\cat ~ A~^/^ —>■ 0, 
beyond which a catastrophic growth of tJ^ begins. 

ii) Two universality classes of passing through the point of singularity. Remarkably, the 
dynamics of passing through the point of singularity can be described analytically (which, 
in itself, is a rather rare event [11]) and it appears quite beautiful. I show analytically 
and conflrm numerically that at the singularity point T = the rate of explosion Qhs 
reaches a maximum fl^g, the dynamics of passing through which proceeds in one of the 
two characteristic scahng regimes, depending on the value of parameter /C = p^/^A/Ac. 

1. Annihilation catastrophe at "frozen" flux (K, ^ 1). In the limit /C ^ 1 at times 
|T| <^Tl = p the particle number n = n^{l + ujq\T\/p + ■ ■ ■) is, in fact, "frozen" n ^ n^, 

and so in the process of explosion in the contracting as Sxl ^ ^\T\/p nearsurface layer 
there arises a sharp gradient of L particles. I show that in this limit up to T -C p the 
flux and the particle number are retained "frozen", hglg — J* and, hence, the explosion 
develops synchronouosly, Qhs = ^Ls = ^s, both before and after the singularity point. 
Equalizing the rates of change of H and L particles diffusion fluxes leads to the scaling 

= 5(r/r;), hjhf^if/h^w{r/Tfi (8) 

where the scahng functions S{C,) = l/y/l + and w{() = C + \/l + C^) the characteristic 
time 7} = 1/^2^, and the amplitudes O.^ , h^ and at the explosion maximum 

= (2//.)Qf = p-'/'^c, if = p"'4aK^. 
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Two striking features of this result arc symmetrical universalization of = 1/|T| beyond 
the scope of interval [—Tf, Tf], 7} ^ and remarkable symmetry T ^ — T, w ^ 1/w. 

Singular "flash" of ilg leads to a singular jump in rj^, the dynamics of which is 
described by the scaling 

[rj'] = [rJ%W{T/Tf), [rj% = c^in^f'/A oc ^l'/'p-'/'A'/^ (9) 

where in accord with the scahng function W{() — C2 /q^ d9/ \/6 [1 + — 9)'^]^/'^ the value of 
[tJ^] reaches at T ~ 7j a maximum max[rj'^] ~ [tJ^]* oc A^^^ — > oo, and then it decays 
as ~ 1 / \/pT. Essentially, that max[T7^] oc p~^^^, so the less is p the more brightly the 
effect is displayed. 

2. Annihilation catastrophe at uniform L particles distribution. Flux breaking effect 
(K, <^1). If A — > oo but p — > so that /C <C 1, the distribution of L species is retained 
uniform {Is = N,Qls = hs) and the final stage of explosion is changed radically. I 
demonstrate that in this case an intermediate scaling takes place 

nHs^^HsQ{T/T*), K^hfvH{T/%), U^lfvL{T/%) (10) 

where the characteristic time % = I/^Ihs amplitudes fijjs^ ^^'^ the explosion 
maximum 

hr = Q^Hs = bHA'/\ If = 6^AV3, 

with q — 2.24, bn — 1-63 and bi — 0.73. It is seen that in the hmit JC <^ 1 the singular 
"flash" of D,Hs occurs much sharper, leading to more abrupt jump in tJ^ 

[rj'] = [rJ%V{T/Z), [rJ% = {q-l)Q^,^A'/^ (11) 

which, in turn, gives rise to a qualitatively new phenomenon: the finite jump of the flux 
by a factor of p:^ 2 for the time \T\break ~ % oc A^^/'^ (flux breaking effect). 

As in the limit /C ^ 1 at T ^ p the relaxation rate rj^ reaches its L-diffusion- 
controUed limit tJ^ ~ ^^o/p, and in the limit /C <S 1 it always grows, I point out finally 
(and demonstrate numerically) that independently of /C in the limit p — > 0, A — > oo 
there arises a most spectacular consequence of the annihilation catastrophe: an abrupt, 
practically instantaneous (on the scale of ujq) disappearance of the flux. 

In conclusion, I believe that annihilation catastrophe may pretend to be not only one 
of the most dramatic manifestations of the reaction- diffusion interplay, but it is a new 
type of catastrophe wherein two explosive processes {^hs and Qls) are developing simul- 
taneously, effectively "compensating" one another so that for an external observer (J) the 
explosion dynamics goes unnoticed up to the point of singularity, in the vicinity of which 
decompensation of explosions is manifested as a sudden singular jump. 
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